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The constraint force equation (CFE) methodology provides a framework for modeling 
constraint forces and moments acting at joints that connect multiple vehicles. With 
implementation in Program to Optimize Simulated Trajectories II (POST2), the CFE 
provides a capability to simulate end-to-end trajectories of launch vehicles, including stage 
separation. In this paper, the CFE/POST2 methodology is applied to the Shuttle-SRB 
separation problem as a test and validation case. The CFE/POST2 results are compared 
with STS-1 flight test data. 


Nomenclature 


A = axial location of the SRB reference point in SSME plume coordinate system 

a^ay,^ = acceleration components along body axes (excluding components due to gravity) ft/sec 2 

Ca - axial force coefficient 

AC a — incremental (proximity) axial force coefficient 

Ci = rolling moment coefficient 

ACi = incremental (proximity) rolling moment coefficient 
C m = pitching moment coefficient 

AC m = incremental (proximity) pitching moment coefficient 
Cv = normal force coefficient 

AC n = incremental (proximity) normal force coefficient 
C„ = yawing moment coefficient 

AC„ = incremental (proximity) yawing moment coefficient 

Cy = side force coefficient 

AC y = incremental (proximity) side force coefficient 

F 'yy F z = aerodynamic forces in axial, lateral and normal directions, lb 
F^con) , f<con) _ j Q | nt constraint force vector for body A, and B 
F^ext) p^ext) _ E xtema i f orce vector acting on body A, and B 
h = altitude, ft 


Ixv, lyzdzx 


L 

M 

MR 

N 


variable index 

moment of inertia about body x, y, and z axes, slugs ft 2 

products of inertia about body x,y, and z axes 

reference length, ft 

rolling moment, lb -ft 

pitching moment, lb -ft 

momentum ratio of BSM jets to freestream 

yawing moment, lb -ft 
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p,q,r = angular velocity components in body axes, deg/s 

q x = freestream dynamic pressure, lb/ft 2 

R = radial location of SRB reference point in SSME plume coordinate system, ft 

S re f = reference area, ft 2 

j(con)^j(con) _ JqJhj. constraint torque vector for body A and B 


V x = freestream velocity, ft/s 

X b ,Y b ,Z b = coordinates in SRB body axes system, ft 

Ax, Ay, A z = relative axial, lateral and normal distances during separation, ft 

a = angle of attack, deg 

[3 = angle of sideslip, deg 

Aa = a(SRB) - a(OET), SRB and OET relative difference in angle of attack, deg 

Ap = p(SRB) - p(OET), SRB and OET relative difference in sideslip, deg 

y = flight path angle, deg 

0,(|),t 1) = Euler angles in pitch, roll and yaw, deg 

p x = freestream density, slugs/ft 3 

a = standard deviation 


Suffixes 

b 

f 

P 


body axes system 

isolated or freestream conditions 

SSME plume impingement 


Acronyms 

ADAMS® = 
AEDC = 
BSM 
CG 

Con Sep = 
ET 

LSRB 

MRC 

NASA = 
NGLT = 
OET 

POST 2 = 

RSRB 
SEE 
SRB 

SSME = 

SSTO 

STS 

SVDS 

TSTO 


Automatic Dynamic Analysis of Mechanical Systems 

Arnold Engineering Development Center 

Booster Separation Motor 

Center of gravity 

Conceptual Separation 

External Tank 

Left Solid Rocket Booster 

Moment reference center 

National Aeronautics and Space Administration 

Next Generation Launch Technology 

Orbiter + ET 

Program to optimize simulated trajectories IT 
Right Solid Rocket Booster 
Synergistic Engineering Environment 
Solid Rocket Booster 
Space Shuttle Main Engine 
Single-Stage-T o-Orbit 
Space Transportation System 
Shuttle Vehicle Dynamic Simulation 
Two-Stage-To-Orbit 


I. Introduction 

The problem of dynamic separation of multiple bodies within the atmosphere is complex and challenging. One 
problem that has received significant attention in literature is that of store separation from aircraft 1 . A similar 
example is the separation of the X-15 research vehicle from the B-52 carrier aircraft. 2 In both of these cases, the 
store and the X-15 vehicle are much smaller in size than the parent vehicle. The other class of stage separation 
problem involves separation of two vehicles of comparable sizes (as in the case of multi-stage reusable launch 
vehicles) where the integrity of each stage is important after separation. A method for computing 6 DOF trajectories 



of hinged/linked vehicles is presented in Ref. 3 and Ref. 4 presents a summary of stage separation capabilities in the 
1960s and 70s. 

NASA studies on stage separation of multi-stage reusable launch vehicles date back to the 1960s. 5 ' 9 More 
recently, Naftel et al 10 " 12 considered staging of two wing-body vehicles. NASA’s Next Generation Launch 
Technology (NGLT) Program identified stage separation as one of the critical technologies needed for the successful 
development and operation of next generation multi-stage reusable launch vehicles 13 . In response to this directive, 
the stage separation analysis and simulation tool called ConSep (Conceptual Separation) was developed as a part of 
this activity. 14 ConSep is a MATLAB® based wrapper to the ADAMS® solver which is an industry standard 
commercial software for solving dynamics problems involving multiple bodies connected by joints. Applications of 
ConSep for two-body (Bimese TSTO) and three-body (Space Shuttle solid rocket booster) separations are discussed 
in Refs. 14 and 15. For performing seamless end-to-end simulations of complete launch vehicle trajectories 
including stage separation, as recommended by the X-43 A (Hyper-X) Retum-To-Flight Review Board 16 , ADAMS® 
has to be linked to a trajectory simulation code like the Program to Optimize Simulated Trajectories (POST) to 
model internal forces and moments at joints that connect multiple vehicles. However, linking ADAMS® (Ref. 17) 
with POST is not advantageous because ConSep uses only a very small fraction of ADAMS® capabilities and it 
significantly increases POST run times for quick conceptual level studies. 

In Refs. 18 and 19, the Constraint Force Equation methodology (CFE) for modeling internal constraint forces 
and moments due to simple joints connecting multiple bodies was developed and its implementation in POST2 was 
discussed. The CFE methodology in POST2 provides a framework for computing the internal constraint forces and 
moments exerted across various joints and applying them as external loads on each vehicle. Thus, the problem of 
motion of multiple bodies connected with joints is reduced to one of multiple free body motion with these additional 
external forces and moments. Consequently, such motions can be simulated with POST2. Thus, with CFE 
implementation, POST2 can be used to perform an efficient end-to-end simulation of launch vehicle trajectories, 
including stage separations. The CFE methodology was tested extensively in a stand-alone mode, using several test 
cases involving fixed, revolute and universal joints. These results were in good agreement with those using 
ADAMS® and Autolev 20 (another industry standard software for simulation of motion of multi-bodies connected 
with joints). The CFE/POST2 implementation was also applied to the separation of the Hyper-X Research Vehicle 
(X^13A) from the Pegasus Booster and the results compared well with the flight test data. 19 

In this paper, CFE/POST2 is applied to a more complex problem of three-body separation, namely that of Space 
Shuttle solid rocket booster (SRB) separation from the Orbiter and External Tank (OET), and the results are 
compared with STS-1 flight test data. 21 ' 23 In this study, only the separation of the SRB from the OET was simulated. 
The end-to-end Space Shuttle trajectory from lift-off to orbit insertion and Orbiter re-entry were not addressed. 

II. Constraint Force Equation Methodology 

To illustrate the basic concept of the CFE methodology, consider the motion of two rigid bodies A and B 
connected by a single joint as shown in Fig. 1. The external forces (Fj^and F^ EXT) ) and torques (T^ ext) and 
Tf-> ) are the resultants of gravity, aerodynamic, and propulsive forces on each vehicle and are routinely 
computed by POST2. Figure 1(b) shows the internal constraint forces and moments (F (CON) and that are 

computed by CFE and act on each vehicle at the joint location. These internal forces and moments constrain the way 
in which one vehicle can move relative to the others, and are dependent upon the external forces acting on each 
vehicle as well as the type of joint. For joints that permit only relative rotation between two bodies, the forces and 
moments on one body are equal and opposite to those acting on the other body, as shown. Finally, the free body 
diagram shown in Fig. 1(c) and the flow chart shown in Fig.2 illustrate the way in which CFE is implemented in 
POST2. At each integration time step, the current version of POST2 computes all the external forces and moments 
acting on each vehicle. The CFE routine works in parallel and computes the internal constraint forces and moments 
required to satisfy the constraints imposed by the joint. Next, these internal joint loads are applied to each vehicle as 
additional external forces and moments, and the POST2 solution is propagated in the usual manner to the next time 
step. Thus, the net external forces and moments on each vehicle are the sum of the usual external forces and 
moments and the joint loads, applied to each vehicle as additional external forces and moments. Consequently, the 
CFE simply augments the vehicle external loads and does not require modification to the POST2 equations of 
motion. The relative displacement between the joint locations is computed at every integration step and is used as a 
measure of how well CFE satisfies the constraint. Additional information on CFE methodology is available in Refs 
18 and 19. 



III. Space Shuttle: Vehicle Description 

The Space Shuttle is designed to launch a variety of payloads into (and/or retrieve from) Earth orbit. The 
Shuttle can deliver up to 65,000 lb for an easterly launch from the Eastern Test Range at Kennedy Space Center 
(KSC) and up to 32,000 lb for a polar orbit launch from the Western Test Range at Vandenberg Air Force Base. The 
Shuttle launch configuration consists of the Orbiter, External Tank (ET) and the two SRBs. A schematic 
arrangement of the launch configuration is shown in Fig. 3. 

The bum time of the SRBs is approximately two minutes. At bum out, the left and right SRBs separate 
simultaneously. The nominal staging of SRBs occurs from 150,000 ft to 180,000 ft altitude and Mach 3.75 to 4.0. 
The SRB lateral and vertical separation is aided by a cluster of four booster separation motors (BSMs) located in the 
forward pods and another cluster of four BSMs on the aft skirts of each SRB. The BSMs provide lateral and normal 
acceleration of the SRBs away from the OET. The relative axial acceleration is achieved by the thrust from the 
SSMEs. After separation, the SRBs are decelerated by parachutes during their descent for splashdown in the ocean, 
approximately 141 nautical miles from the launch site. The SRBs are recovered and are refurbished for reuse. A 
schematic illustration of SRB separation, descent and recovery is presented in Fig. 4. 

The SRBs are attached to the ET at one point on the forward end and at three points on the aft end. The forward 
attachment point is used to transfer the thrust loads from the SRB to the ET. The aft attachment points consist of 
three stmts. At SRB separation, the bolts holding the forward and aft attachment points are severed, the aft joints are 
released, and the SRBs rotate about the forward joints by approximately 1 deg prior to release. Simultaneously, the 
forward and aft BSMs on each SRB are fired. The BSMs bum for about 0.7 sec. The SRB separation is complex and 
is influenced by several factors. The OET cannot make a lateral maneuver to aid separation because an SRB is 
coming off each side. Consequently, the OET is flown straight and in attitude-hold mode. Because the Orbiter 
engines are significantly above the OET center of gravity, they need to be canted for achieving moment balance for 
attitude-hold mode as shown in Fig. 5(a). As a result, the OET moves in z-direction towards the SRBs. The BSMs 
have sufficiently high thrust to move the SRBs away from the OET, resulting in a safe separation. 

The nominal mass and inertia at SRB separation for STS-1 are given in Table I (Refs. 21-23). The center of 
gravity locations and aerodynamic moment reference points are shown in Fig. 5. During SRB separation, the thrust 
produced by each SSME was assumed to be equal to 467,000 lb, and the net thrust was assumed to pass through the 
center of gravity of the OET to reflect the “attitude hold” mode of the OET. The canting angle of the net thrust was 
estimated to be 13.35 deg as shown in Fig. 5(a). 

The flight parameters at SRB separation for STS-1 are presented in Table II. The SRB propulsion characteristics 
for STS-1 are given in Table III. It may be noted that Ref. 23 gives the thrust values of 18,297.0 lb (left SRB) and 
24049.0 (right SRB) for STS-1 at separation but time history of thrust tail-off is not presented. On the other hand, 
Ref. 23 gives the time history of nominal SRB thrust tail-off for the SRB separation sequence and these data are 
used in this paper. The time history of nominal thrust tail off from the initiation of SRB separation (with a 6 sec 
delay from separation cue) is shown in Fig. 6. The BSM nominal locations are presented in Table IV. Each of these 
BSM jets are nominally inclined 40 deg in the pitch plane and 20 deg in the roll-yaw plane as shown in Fig. 7. The 
variations in BSM cluster orientations from their nominal values and the cluster propulsion characteristics for STS-1 
are given in Table V. 


IV. Aerodynamic Database 

The Shuttle SRB separation aerodynamic database is comprised of three components: (i) an isolated 
aerodynamic database for each vehicle when they are in freestream conditions, (ii) proximity coefficients 
(increments in relation to freestream coefficients and including mated conditions) for each vehicle, and (iii) plume 
impingement forces and moments on the SRBs to account for SSME plume impingement when that happens. For 
this analysis, the isolated and plume impingement aerodynamic characteristics for the left and right SRBs were 
assumed to be identical with proper sign adjustment. Note that, when the SRBs are not subject to SSME plume 
impingement, these forces and moments are zero. 

The relative orientations and location during separation are defined in Fig. 8. Here, Ax, Ay, and Az are 
orthogonal displacements of SRB from its mated position. The variable Ax is measured positive aft, Ay is positive 
outboards, and Az is positive downward. For the mated position, Aa = A(3 = Ax = Ay = Az = 0. The relative angles of 
attack and sideslip are defined as, Aa = a(SRB) - a(OET) and A(3 = p(SRB) - P(OET). 

The proximity aerodynamic coefficients depend on eight variables a, (3, Aa, A(3, Ax, Ay, Az and MR. The BSM 
jet momentum ratio (MR) was evaluated using the data given Table V and Figs. 8 and 9. The MR was found to be in 



excess of 730. However, the hypercube proximity aero data goes up only to MR = 222.3. Therefore, when 
MR>222.3, using the “no extrapolation” approach recommended in Ref. 21, data for MR = 222.3 was used. At 
BSMs bum out, MR=0 and the power-off proximity data was used. Since the staging maneuvers last only for few 
seconds, the Mach number is usually assumed constant. The SSME plume impingement forces and moments depend 
on four variables A, R, 0, and op as shown in Fig. 10. Here, A and R are the axial and radial locations of the SRB 
reference point in SSME plume coordinate system, and 0 and are the pitch and yaw angles. The geometrical shape 
of the SSME engine plume was modeled using the data on the jet boundary given in Ref. 23. The plume 
impingement data 19 covers axial distances only up to 1700 in. For axial distances greater than 1700 in, the last value 
was held constant. The development of a proximity aerodynamic database for the Shuttle SRB separation was a 
difficult and complex task, because the number of data points needed to populate a conventional database with eight 
independent variables would be unusually large. To get around this difficulty, the Shuttle program developed a novel 
method called “hypercube 24 , which used far fewer data points. This method allowed data to be placed only along 
required separation paths. For each selected Ax, two 4-dimensional hypercubes were situated to cover all anticipated 
trajectory values in Ay, Az, Act, and A (3. The inner hypercube was designed to include nominal separation trajectories 
with 3a dispersions, and the outer cube was designed to encompass all dispersions, including SSME failures. These 
hypercubes were not constrained to have parallel opposite sides. A nominal trajectory through hypercube space is 
schematically illustrated in Fig. 11. The range of variables in hypercube data is shown in Table VI and aerodynamic 
reference parameters are listed in Table VII. Additional information on the hypercube concept may be found in 
Ref. 19. The aerodynamic forces and moments during separation are defined as 


V^C^+AC^+F^ 

(1) 

V’S„ f w(c 0M + AC 0 ,4 + T pai 

(2) 


where, C Qt4 and AC Qt4 are any one of the three isolated or proximity force coefficients (axial, normal or side force) 
C Qt4 and AC Qt4 are any of the three moment coefficients(pitching, rolling and yawing) and, i = 1,2,3 for right SRB, 
OET and left SRB respectively. 

To verify proper implementation of the hypercube interpolator, Ref. 24 provides one check case. The input 
variables for this check case are: a = 1.0 deg, (3 = 4.0 deg, Ax = 186.0 in, Ay = 65.0 in, Az = 133.0 in, Act = -4.4 deg, 
A(3 = -0.9 deg, and MR = 196.1634. Present results are in good agreement with that check case as shown in Table 
VTII. This exercise indicates that both the hypercube interpolator routine and the proximity aerodynamic coefficients 
were properly implemented in this study. However, it may be noted Ref. 24 specifies MR = 200 for this check case 
but a detail verification indicated that the value of MR=196.1634 better matches all the other parameters for this 
check case including the aerodynamic coefficients. 

The ConSep/ADAMS results (Ref 15) were based on modeling the proximity aerodynamic coefficients, using a 
response surface method, because the hypercube interpolator was not available at that time. Since a hypercube 
interpolator is the recommended approach, the Con Sep results included in this paper were recomputed using the 
hypercube interpolator. 

The isolated aerodynamic coefficients for each vehicle and the plume impingement forces and moments were 
obtained from the Space Shuttle Launch vehicle databook (Ref. 24). The proximity data tables and the hypercube 
interpolator routines were provided by Hoff and Paget (Ref. 25). 

The CFE/POST2 simulation inputs include the following: (i) definition of body-fixed 
coordinate systems for each vehicle in relation to the selected ground-fixed (inertial) system, (ii) 
mass, inertia, and center of gravity of each of the three vehicles in respective body-fixed 
coordinate system, (iii) initial altitude, velocity and flight path angle, (iv) initial angle of attack 
of each vehicle, (v) proximity aerodynamic data tables and hypercube interpolator routines, (vi) 
isolated (freestream) aerodynamic coefficients, (vii) plume impingement forces and moments, 
(viii) aerodynamic reference parameters for each vehicle, (ix) thrust magnitudes and orientations 
for each of the forward and aft BSMs for both right and left SRBs, (x) thrust magnitudes for each 
SSME. The time/event for the release of front and aft joints, and the time for starting and 



stopping the integration of equations of motion and step size, are presented in Table IX. Each aft 
joint was modeled as a fixed joint, released instantaneously, and each forward joint was modeled 
as a revolute joint that released when the SRBs rotated by 1 deg with respect to mated condition. 
The effect of forces in the attachment hardware at SRB release were not included in CFE/POST 
2 simulation because the current version of CEE does not support preloaded joints. 

The results of CFE/POST for STS-1 SRB separation are presented in Figs. 12 to 18 along with STS-1 flight data, 
Rockwell/Boeing SVDS (Shuttle Vehicle Dynamic Simulation), and ConSep simulation results. 15,22 The predicted 
relative locations of the vehicles during separation are in good agreement with the flight data (Fig. 12). However, 
differences exist in the lateral and vertical separation distances Ay and Az. In the present ConSep and CFE/POST2 
simulations, the estimated SSME thrust orientations (Fig. 5a) were used and it is possible that for STS-1 flight, there 
could be variations from these values and these factors directly contribute to observed differences in Ay and Az. The 
component accelerations for OET and for the left and right SRBs are in good agreement with the flight data (Figs. 
13 to 15). However, the component angular velocities show differences from the flight data (Figs. 16 to 18) and the 
reasons for these differences are not clear. The gimbal angles of each of the three SSME were continuously varied 
during separation in STS-1, and these variations were not modeled in CFE/POST and ConSep. Furthermore, strut 
forces acting at the aft joint at release were not included, and, the rolling moment coefficient for each SRB was set 
equal to zero. It is possible that such factors contribute to these differences. 

The Synergistic Engineering Environment 26 (SEE) was used to create animations of the staging maneuvers. The 
SEE used the geometry models of the SRBs, OET, and the time history of trajectory variable to generate these 
animations. The BSM plume shapes were computed using the data of Ref. 21 and were included in SEE animation 
for the duration of their bum (0.68 sec). The SSME plumes are on for the entire 4 sec duration. Some sample 
snapshots of SEE animation are presented in Fig. 19. It is interesting to note that the SSME plume impinges on the 
SRBs for t > 3.0 sec. 


V. Concluding Remarks 

The constraint force equation (CFE) methodology provides a framework for modeling constraint forces and 
moments acting at joints that connect multiple vehicles. With implementation in POST II, the CFE provides a 
capability to simulate end-to-end trajectories of launch vehicles including stage separation. The POST/CFE 
simulation methodology was applied to the STS-1 SRB separation as a further test and validation case for 
CFE/POST2. The results compare fairly well with flight data and this exercise serves as a test/validation for the 
CFE/POST 2 algorithm. With this validation, the CFE/POST 2 can be used for end-to-end simulation of launch 
vehicle trajectories, including stage separation, for system level studies. 
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Table I. Vehicle Mass and Inertia Properties at Staging (Ref. 22) 


Configuration 

Weight, 

lb 

Ixx, 

Slugs ft 2 

iyy, 

Slugs ft 2 

Izz, 

Slugs ft 2 

Ixy, 

Slugs ft 2 

Iyz, 

Slugs ft 2 

Ixz, 

Slugs ft 2 

Orbiter/+ET 

1,473,332 

5,241,320 

83,826,248 

80,221,300 

124,263.3 

-12,024,025.0 

35,820.9 

Left SRB 

180,881 

192,915 

11,348,219 

1 1 ,35 1 ,222 

-27,923.0 

7,508.0 

688.0 

Right SRB 

181,833 

241,702 

11,463,621 

11,466,376 

24,651.0 

9,670.0 

-1,043.0 


Table II. Flight Parameters at SRB Separation for STS-1 (Ref. 22) 


Parameter 

Value 

Altitude, ft 

173,857 

Velocity, ft/sec 

4112.0 

Mach Number 

3.877 

Dynamic Pressure, lb/ft 2 

12.83 

Angle of Attack, deg 

3.0 

Sideslip, deg 

0.75 

Flight Path Angle, deg 

36.966 

Body Roll Rate, deg/sec 

-0.3 

Body Pitch Rate, deg/sec 

0.4 

Body Yaw Rate, deg/sec 

0 


Table III. SRB Propulsion Characteristics at Separation, STS-1 (Ref. 22) 



Parameter 

Left SRB 

Right SRB 



Propellant Mean Bulk 

as n 

68U) 



Temperature (PMBT), deg F 




Pitch gimbal angle, deg 

0.042 

0.064 



Yaw gimbal angle, deg 

-0.679 

0.672 



Table IV. BSM Nominal Locations in SRB coordinates (Ref. 22) 

BSM No 

Type 

Xb (inches) 

Yb (inches) 

Zb (inches) 

l 

Forward 

290.230 

19.625 

26.616 

2 

Forward 

303.888 

20.811 

30.863 

3 

Forward 

290.23 

2.351 

32.772 

4 

Forward 

303.888 

3.896 

37.019 

5 

Aft 

1886.722 

59.755 

82.245 

6 

Aft 

1886.722 

40.573 

93.418 

7 

Aft 

1886.722 

28.357 

97.864 

8 

Aft 

1886.722 

15.548 

100.682 


Table V. Variations in BSM Orientations and Propulsion Characteristics at Separation (STS-1) 


Parameter 


Left SRB 
Forward 


Left SRB 
Aft 


Right SRB 
Forward 


Right SRB 
Aft 


A(j), deg 


0.04 


0.113 


0.18 


0.05 











A0, deg 


0 


0.10 


0.14 


0.071 


Maximum Chamber 
Pressure, psia 

1844.9 

1832.1 

1854.7 

1848.4 

Maximum Thrust, lb 

92,096.0 

91,457.0 

92,587.0 

92,272.0 

Average Thrust A lb 

89,620 

88,162 

90,776 

90,884 

Average Duration, sec 

0.692 

0.696 

0.671 

0.686 


Table VI. Range of Independent Variables in Hypercube Proximity Aerodynamic Database 

Plume-off (MR = 0) 


a (deg) 

P (deg) 

Ax (inches) 

Ay (inches) 

Az (inches) 

Aa (deg) 

A(3 (deg) 

-10 to 10 

-10 to 10 

0 

0 

0 

0 

0 

-10 to 10 

-10 to 10 

100 

10 to 110 

40 to 250 

-7.0 to 0 

-5.5 to 1.0 

-10 to 10 

-10 to 10 

300 

30 to 260 

40 to 550 

-17.0 to 0 

-15.0 to 2.0 

-10 to 10 

-10 to 10 

600 

90 to 510 

140 to 800 

-30.0 to 0 

-20.0 to 3.0 

-10 to 10 

-10 to 10 

1100 

100 to 700 

180 to 900 

-33.0 to 0 

-20.0 to 3.0 

-10 to 10 

-10 to 10 

1700 

200 to 800 

300 to 1000 

-34.0 to -5.0 

-20.0 to 8.0 

Plume-on (MR = 1334, 2223) 

a (deg) 

P (deg) 

Ax (inches) 

Ay (inches) 

Az (inches) 

A a (deg) 

A(3 (deg) 

-10 to 10 

-10 to 10 

0 

0 

0 

0 

0 

-1 0 to 10 

-10 to 10 

100 

10 to 110 

40 to 250 

-7 to 0 

-5.5 to 1.0 

-10 to 10 

-10 to 10 

200 

20 to 150 

60 to 280 

-7 to 0 

-6.5 to 0.5 








Table VII. Aerodynamic Reference Parameters 


Vehicle 

Reference Area 

Reference Length 

Moment Reference Point 

Right SRB 

2690.0 sq ft 

107.525 ft 

X=88.27, Y=0. Z=0 in SRB 
coordinate system with origin at 
nose 

OET 

2690.0 sq ft 

107.525 ft 

X=64.54 ft, Y=0, Z= 4.17 ft, in 
OET coordinate system with 
origin at nose 

Left SRB 

2690.0 sq ft 

107.525 ft 

X=88.27, Y=0. Z=0 in SRB 
coordinate system with origin at 
nose 


Table VIII. Verification of Hypercube Interpolator Implementation 


Aerodynamic 

Coefficient 

Present Result 

Check Case (Ref. 24) 

CNoet 

-0.08119 

-0.08117 

CmoET 

-0.07141 

-0.0714 

CYoet 

0.03983 

0.03984 

CnoET 

-0.03674 

-0.03674 

CloET 

0.01277 

0.01277 

CNsrb 

0.00399 

0.00399 

Cm S RB 

0.00899 

0.00899 

CYsrb 

0.011 

0.011 

Cn S RB 

-0.0009 

-0.0009 


Table IX. CFE/POST 2 Simulation Parameters 


Parameter 

Value 

Simulation time, sec 

4.0 

Integration step size, sec 

0.001 

Aft (fixed) joint release, sec 

0.02 

Front (re volute) joint 
release, deg 

1.0 




















a) External forces. b) Internal forces. c) Free body diagram. 

Figure 1. Decomposition of external forces and internal constraint forces to free body diagram solved by CFE 

method. 








Figure 3. Space Shuttle vehicle configuration. (Ref. 20). 





Figure 4. Schematic illustration of Space Shuttle staging events. 
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(a) Center of gravity locations 
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(b) Aerodynamic moment reference points. 


Figure 5. Center of gravity and aerodynamic moment reference points for STS-1. 
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Figure 7. BSM cluster locations and orientations. 





Figure 8. Relative orientations and locations during separation. 
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Figure 9. Variation of BSM jet momentum ratio with freestream dynamic pressure (Ref. 20). 
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Figure 10. Schematic illustration SSME plume impingement coordinate system (Ref. 24). 



Figure 11. Hypercube data structure and nominal trajectory (Ref. 24) 



Separation Distances 




r — cr ' '4 cr " 


ConSep LSRB 


ConSep RSRB 


O Flight LSRB 


E Flight RSRB 
* SVDS LSRB 


+ SVDS RSRB 

4 

— POST2 CFE LSRB 


POST2 CFE RSRB 




Figure 12. Variation of relative distances Ax, Ay and Az with time during SRB separation, STS^l. 
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Figure 13. OET: Variation of acceleration components during separation, STS-1. 
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Figure 14. Left SRB: Variation of acceleration components during separation, STS-1 
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Right SRB Accelerations (excluding gravity) 
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Figure 15. Right SRB: Variation of acceleration components during separation, STS-1 
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Figure 16. OET: Variation of angular velocity components during separation, STS-1 
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Figure 17. Left SRB: Angular velocity components during separation, STS-1 
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Figure 18. Right SRB: Angular velocity components during separation, STS-1 
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(a) t = 0. 
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(b) t=0.5 sec 

Sim: Time = 1.0 Seconds 
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(d) t=1.5 sec 


Sim: Time = 2.0 Seconds 



(e) t = 2.0 sec. 


Sim: Time ■ 2.5 Seconds 


I. BZS src 


Figure 19. Snapshots of SEE animation 



